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1.1 INTRODUCTION 

Harmonic generation is one of the earliest discovered and studied nonlinear opti- 
cal processes. For 40 years, since the first experimental demonstration of second- 
harmonic generation (SHG) by Franken and co-workers [1] followed by its rigorous 
theoretical description by Bloembergen and Pershan [2], the harmonic generation has 
unceasingly been attracting much attention [3]. In particular, harmonic generation 
has been apphed as a source of nonclassical radiation (see references [4, 5] for a de- 
tailed account and bibliography). It was demonstrated that photon antibunched and 
sub-Poissonian light [6, 7], as well as second [8] and higher order [9, 10] squeezed 
hght can be produced in SHG. In experimental schemes, second-harmonic generation 
is usually applied for the sub-Poissonian and photon-antibunched light production, 
whereas second-subharmonic generation (also referred to as the two-photon down 
conversion) is used for the squeezed-light generation [4, 1 IJ. Non-classical effects 

tXo appear in Modern Nonlinear Optics, ed. M. Evans, Advances in Chemical Physics, vol. 1 19(1) (Wiley, 
New York, 2001). This is a part of the chapter on Nonlinear phenomena in quantum optics by J. Bajer, 
M. Dusek, J. Fiurasek, Z. Hradil, A. Luks, V. Pefinova, J. Rehacek, J. Pefina, O. Haderka, M. Hendrych, 
J. Pefina, Jr., N. Imoto, M. Koashi, and A. Miranowicz. 
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in higher-harmonic generation have also been investigated, including sub-Poissonian 
photocount statistics [5, 7, 12, 13], squeezing [5, 14, 15], higher-order squeez- 
ing [16, 17] according to the Hong-Mandel definition [9] or higher-power-amplitude 
squeezing [18, 17] based on Hillery's concept [10]. In this contribution, we will 
study photocount statistics of second and higher harmonic generations with coherent 
light inputs. 

Photocount noise of the observed statistics can simply be described by the (quan- 
tum) Fano factor [22] 



(1.1) 



where (n) is the (ensemble) mean number of detected photons and ((An)^) is the 
variance of photon number. We also analyze the global (quantum) Fano factor 
defined to be [23]: 

- ((n)) ((n)) ' 

where the mean values ((n*^)) are obtained by the ensemble and time averaging, i.e., 

{{h'))= Mm i f{h\t))dt. (1.3) 
In classical trajectory approach, the Fano factor is defined to be 



pS^iAnl^r^^ (1.4) 

n n 

as a semiclassical analogue of the quantum Fano factor. The mean values n'' in Eq. 
(1.4) are obtained by averaging over all classical trajectories as wiU be discussed in 
detail in Sects. 1.2.3 and 1.3.3. 

Coherent (ideal laser) light has Poissonian photon-number distribution thus de- 
scribed by the unit Fano factor. For F < 1, the light is referred to as sub-Poissonian 
since its photocount noise is smaller than that of coherent Ught with the same inten- 
sity. Whereas for F > 1, the light is called super-Poissonian with the photocount 
noise higher than that for coherent light. 

We shall compare different descriptions of photon-number statistics in harmonic 
generation within quantum, classical and semiclassical approaches. First, we will 
study the exact quantum evolution of the harmonic generation process by applying 
numerical methods including those of Hamiltonian diagonaUzation and global char- 
acteristics. As a brief introduction, we will show explicitly that harmonic generation 
can indeed serve as a source of nonclassical light. Then, we will demonstrate that the 
quasi- stationary sub-Poissonian Ught can be generated in these quantum processes 
under conditions corresponding to the so-called no-energy-transfer regime known 
in classical nonlinear optics. By applying method of classical trajectories, we will 
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demonstrate that the analytical predictions of the Fano factors are in good agreement 
with the quantum results. On comparing second [19], third [20] and higher [21] 
harmonic generations in the no-energy-transfer regime, we will show that the highest 
noise reduction is achieved in third-harmonic generation with the Fano-factor of the 
third harmonic equal to « = 13/16. 

1.2 SECOND-HARMONIC GENERATION 
1.2.1 Quantum analysis 

The quantum process of second-harmonic generation (SHG) can be described by the 
following interaction Hamiltonian [4, 5]: 

H = hg (afa^ + a\^a2^ , (1.5) 

where oi and 02 denote annihilation operators of the fundamental and second- 
harmonic modes, respectively; 5 is a nonlinear coupling parameter. The Hamiltonian 
(1.5) describes a process of absorption of two photons at frequency loi and simul- 
taneous creation of a new photon at the harmonic frequency uj2 = 2oji, together 
with the inverse process. Unfortunately, no exact solution of quantum dynamics 
of the model, described by (1.5), can be found. Thus, various analytical approxi- 
mations or numerical methods have to be applied in the analysis of the conversion 
efficiency, quantum noise statistics or other characteristics of the process [5]. Due to 
mathematical complexity of the problem, the investigations of nonclassical effects in 
harmonic generation have usually been restricted to the regime of short interactions 
(short optical paths or short times). Theoretical predictions of quantum parameters 
(including the Fano factor or, equivalently, the Mandel Q-parameter) were obtained 
under the short time approximation only (see, e.g., [4, 5, 13]). This is a physi- 
cally sound approximation in case of weak nonlinear coupling of optical fields. The 
Fano factors under the short-time approximation (i.e., for gt <^ 1) for coherent 
inputs ai = ri exp (icpi) and a2 — r2 exp (^^2) are given by the expansions (for 
ri,r2 ^0): 

F^ = l-4siner2ffi 

+ {4rrVi - 2rf + 8[2 + cos(20)]ri} {gtf + 0{{gtf}, 

F^ = l-^sinO rlr^igtf 

+^ {2ri + Urfrl - [4 + 3 cos(20)]r^} {gt)^ + 0{{gtf}, (1.6) 

where 6 = 2([)i — (t>2 and 0{x} denotes the order of magnitude. Eq. (1.6) determines 
whether the generation of harmonics (ui + co ^ 2lo) or subharmonics (2w ^ ui + ui) 
occurs. It also determines the sub-Poissonian or super-Poissonian photon-number 
statistics of Ught generated during the short-time interactions. For spontaneous SHG 
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Fig. 1.1 Fano factors of the fundamental, F^, and the second-harmonic mode, F^, in the 
long-time interaction for initial coherent states with real amplitudes (a) ai = 6,a2 = 1, 
and (b) ai = 6, 0:2 = 3. Case (a) is a typical example of super-Poissonian behavior in both 
modes outside the no-energy-transfer regime. In case (b), the harmonic mode exhibits stable 
sub-Poissonian statistics with ~ 0.83. It is a characteristic example of the sub-Poissonian 
behavior within the no-energy-transfer regime along the line |ai| = 2|a2|- 



process (i.e., for r2 = 0), the well-known expansions for the quantum Fano factors 

are = l - 2{ngtr + -rf {Srj + l) {gtf + 0{{gtf}, 

= 1~ j{ngtf + ^rt (SGrl + 17) igtf + 0{igtf}, (1.7) 
or, equivalently, for the normally-ordered photon-number variances [6, 7, 4]: 
(: (Ani)' :) ^ {{Ah^f) - (ni) = -2rt{gtf + 0{{gt)'}, 
{■.{Amf:) ^ {{Ah2f)-{h2) = -^rt{gtf + 0{{gtf}. (1.8) 

It is seen that the photon-number statistics of fundamental mode exhibits, in the 
short-time regime, much stronger sub-Poissonian behavior than that of harmonic 
mode. 

For longer interaction times (gt > 1), there are no exact analytical solutions, 
thus the numerical analysis has to be applied. We have used two methods to study 
the quantum dynamics: (i) the well-known Hamiltonian diagonalization proposed 
by Walls and Barakat [24] and (ii) the method of global characteristics based on 
manipulation with spectra [23]. These methods can be applied for arbitrary initial 
photon statistics. Nevertheless, for the purpose of our paper, we restrict our analysis 
to the initial coherent fields solely. Due to computational difficulties, the results can 
be obtained for small numbers of interacting photons only. The analysis of about one 
hundred interacting photons reaches practically the computational capabiUties of the 
standard mathematical software. 

Analysis of a typical evolution of the Fano factors F^2' such as presented in Fig. 
1.1(a), leads to the conclusion that after initial short-time (gt < 1) relaxations in 
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Fig. 1.2 Global Fano factor, Fj', of the second-harmonic mode as a function of the initial 

coherent state amplitudes ai and a2 with = 0. It is seen that the harmonic mode exhibits 
globally sub-Poissonian behavior (F^ < 1) near the diagonal |ai| = 2\a2 \ wdO = 0. The 
darker region the higher value of • Th^ counter lines are drawn at = 1,1. 5, •••,5. 5. 



both modes, a strongly super-Poissonian (-Fi 2 ^ !)• This behavior occurs for the 
majority of initial coherent states \ai) and \a2) except a certain set of initial states 
concentrated along the line |ai | = 2\a2\ > and ~ [see Fig. 1.1(b)]. The same 
conclusion can be drawn by analyzing the global Fano factors F^2- We find that the 
global Fano factor of the harmonic mode remains independent of amplitude | and 
equal to = 0.83 < 1 along the line |ai| = 2\a2\ (see Figs. 1.2 and 1.3). As 
depicted in Fig. 1 . 1 (b), when the initial relaxation oscillations fade out, the harmonic 
mode remains sub-Poissonian for a long interaction time interval. In the classical 
theory of SHG, this case is referred to as the no-energy- transfer regime [25], because 
of the conservation of energy in every mode. We have found a quantum analog of this 
regime for coherent inputs with amplitudes satisfying the conditions: |ai| = 2|a2| 
and 0-0. 

By analyzing Figs. 1.2 and 1.3, an intriguing question arises: Why does the 
harmonic-field photocount statistics in the no-energy-transfer regime remain sub- 
Poissonian with the Fano factor almost independent of the interaction time gt > 17 
This behavior can be understood better by plotting the Husimi Q-function. Let us 
have a look at the snapshots of typical evolution of Q-functions for both modes at 
six times gt = 0, 0.5 • • • , 2.5 with initial amplitudes ai = 6, 02 = 3 (Fig. 1.4). 
These results were obtained numerically and represent the exact quantum solution of 
the model (1.5). One can observe how the cross-sections of the Q-functions change 
from circles (for initially coherent fields) through crescents into rings. We note that 
both modes have relatively small photon-number variances and small Fano factors, 
Ri 1.50 and F^ w 0.83 [see also Fig. 1.1(b)]. The ring shapes, once formed. 
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Fig. 1.3 Global Fano factors and along the line |ai | = 2|a2| for ^ = 0. Dotted 
lines denote the RMS deviation of oscillations in the long-time interaction. It's seen that the 
harmonic mode exhibits globally sub-Poissonian behavior. 

are very stable. So, not only the Fano factors, but the entire quantum states become 
stationary. 

The Q-functions are very wide, thus no linearization of the quantum problem is 
possible and no pure quantum technique can be used for estimation of the observed 
values « 1.50 and F^ « 0.83. However, good quantitative explanation of these 
numerical values can be obtained by the method of classical trajectories as will be 
shown in Sect. 1.2.3. 

Our discussion is focused on photon-number statistics rather than squeezing or 
other phase-related properties. Nevertheless, by analyzing the Q-function evolution 
presented in Fig. 1 .4, we can draw the conclusion that squeezing cannot be observed 
for initial coherent fields at interaction times exceeding the relaxation time. In fact, 
the quadrature squeezing variances (fc = 1, 2) 

^ HAX.r) = ([A(a,e-'^ + 4e'^)]2) (1.9) 

are monotonically rising from the standard shot-noise-limit (S'Q = 1) to much more 
noisy state with the saturated quadrature variances 

S'l^ = 1 + > 1, 

= 1 -h 2r2 > 1. (1.10) 

It is evident that the no-energy-transferregime is not useful for the quadrature squeez- 
ing generation. 

1.2.2 Classical analysis 

Complete quantum solution of the model given by Hamiltonian (1.5) can be found by 
applying sophisticated numeric methods on a fast computer only. However, since we 
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Fig. 1.4 Quantum evolution of the Q-function for the fundamental (outer contour plots) and 
the second-harmonic mode (inner plots) at six time moments for initial coherent states with 
ai — Q, a2 — A, — 0. Solution obtained by quantum numerical method. 



are interested in a special type of solution for strong fields, we can adopt approximate 
classical and semiclassical methods to obtain some analytical results. 

In analogy with Eq. (1.5), the classical model of SHG can be described by 

g[ala*2 + {alfa2], (1.11) 

where ai and a2 are complex amplitudes of the fundamental and second-harmonic 
modes, respectively, and is a nonlinear coupling parameter. The exact solution of the 
model, described by (1.11), is well-known (see, e.g., [26]). The solution is periodic 
and can be written in terms of the Jacobi elliptic function. A few special cases 
(e.g., the phase-matched second-harmonic generation) have monotonous solution 
described by hyperbolic functions. The classical solution is a good approximation 
for strong fields, for which gives correct predictions of the output light intensities 
and frequency-conversion efficiency. Unfortunately, it cannot be used to describe the 
photocount noise and other statistical properties of generated light. Now, we will 
summarize some classical results, which we will be used in the method of classical 
trajectories. 

The Hamiltonian (1.11) for the classical SHG leads to the following system of 
complex differential equations [26] 



di = 



-2\gala2, 



8 



da = -igaf. (1.12) 

One obtains, after substitution of ak = rkS^'^'', a new system of real equations for 
the amplitudes and phases: 

ri = —2rir2smd, 
r2 = r^sinff, 

= (r^/ra -4r2) cos 0, (1.13) 

where 9 = 20i — 02- The system has two integrals of motion: E = r\ + 2r| = 
ni + 2n2 and F = r\r2 cos 6. By extracting n and 6 from Eq. (1.13), we get the 
following equation for r2: 

{r2f2/gf + T^=rl{E-2rl)\ (1.14) 
or even in simpler form for the intensity n2 = r^: 

{n2/2gf = n2{E -2n2f . (1.15) 
Separation of t and n2 leads to the equation 

2gdt = (1.16) 

{E-2n2f-T^ 

which can be rewritten as 

Agdt = —== ^'^^ (1.17) 
^J{a - na) (6 - na) (n2 - c) 

where the numbers a, 6, c are the roots of cubic equation n2 {E — 2n2)^ — = 0. 
For c <u <h < a, the solution of Eq. (1.17) reads as 



r 



dx 2 f u — c 



asn 



ur^M (1.18) 



in terms of the inverse Jacobi elliptic function, asn {x, k), with parameter fc = y ^rf- 
Finally, the inversion of (1.18), gives the required solution 

n2{t) =c+{b-c) sn'^[2g^/^^{t-to),k] , (1.19) 

where sn(u, k) is the Jacobi elliptic function with the same parameter k. Solution 
(1.19) can be simphfied in special cases. In particular, the well-known elementary 
solution is obtained for second-harmonic generation from vacuum, where ri (0) = r 
and r2 (0) = 0. In this case k = 1 and the Jacobi elUptic function simpUfies to 



SECOND-HARMONIC GENERATION 9 



hyperbolic tangent. Thus, the solution reads as 

ri{t) = rcosh{V2rgt), 

r2{t) = -^t&nh{V2rgt) (1.20) 
v2 

and 6{t) = n/2. Subharmonic generation does not occur in this classical model, 
since for ri (0) = and r2 (0) = r implies that ri (t) = 0. r2 (t) = r for any 
evolution time t. Another important special case of solution (1.19) can be obtained 
for the initial zero phase difference, 9 (0) = 0, and the initial amplitudes satisfying 
n (0) = 2r and r2 (0) r. Here, E = Sr^, T = 4r^, a = 4r'^, b = c = r'^, k = 
and Jacobi elhptic function simplifies to trigonometric sinus. Finally, this elementary 
solution reads as 

ai (t) = 2re-2''-9*, 

a2{t) = re-^'''9*, (1.21) 

which corresponds to the no-energy-transfer regime, in which energy is conserved 
in every mode. Phase trajectories of that solution are presented in Fig. 1.5(a). The 
slightly perturbed solution in the no-energy-transfer regime can also be approximated 
by fc « and elementary function sinus with small amplitude (6 — c) ~ [see Fig. 
1.5(b),(c)]. 



1.2.3 Classical trajectory analysis 

The answer to our question concerning the origin of sub-Poissonian behavior can 
be found by the method of classical trajectories. The method is very general. It 
can be applied in the analysis of almost every nonlinear quantum process. Even 
external pumping and energy losses can be easily described. In the classical trajectory 
approach to SHG [27], deterministic solutions of the classical SHG are used, while 
quantum noise of initial fields is artificially simulated by Gaussian distribution. One 
can study the time evolution of the bunch of trajectories Hke the evolution of quantum 
distributions. This semiclassical method can often shed some light on complicated 
quantum dynamics. For strong inputs, where the quantum noise can be assumed 
small, the method gives surprisingly good results. 

According to the classical trajectory method, one assumes that the input stochastic 
amplitudes are of the form 

ai = ri -|-a;i -l-iyi, 

a2 = r2+X2 + iy2, (1.22) 

where rk are coherent complex amplitudes, whereas Xk and yk are real and mutually 
independent Gaussian stochastic quantities with identical variances 



a2 = l/4«r^ 



(1.23) 




flea flea 



M (d) 



s 




Fig. 1.5 Classical trajectories of the fundamental and second-harmonic modes in the no- 
energy-transfer regime for ai = 6 and a2 = 3: (a) classical evolution according to Eq. 
(1.19) for < gt < 0.5; (b) fifty random trajectories out of 10,000 trajectories used in the 
simulations for < gt < 0.5; (c) same as in (b), but for < gt < 5; (d) snapshot of the 
Q-function, obtained from 10,000 random trajectories at time gt = 5.0. Triangles denote 
starting points {gt = 0) of the trajectories and circles are their ends. 



By analogy with our quantum analysis, we calculate the semiclassical Fano factor, 
defined by Eq. (1.4), and quadrature squeezing variance 

Si = {KX^ = [A(afee-'^ + a^e'^)]^ (1.24) 

as counterparts of quantum parameters (1.1) and (1.9), respectively. By applying the 
method of classical trajectories with the noise variance given by Eq. (1.23), we find 
the semiclassical quadrature squeezing and Fano factor given by: 

St = 4<t2 = 1, 

« W ^ 1, (1.25) 

respectively. According to the described method, one needs to solve thousands of 
the classical SHG trajectories. The mean values are simply obtained by averaging 
over all these trajectories. In Fig. 1.6, we have presented graphically snapshots in 
a selected time-interval of all complex solutions in phase space. These clouds of 
points naturally correspond to the Q-functions in the quantum picture (see Fig. 1 .4). 
We have found that this semiclassical method gives the results surprisingly similar 
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to the quantum results even for relatively weak fields! This very good agreement is 
clearly seen by comparing Figs. 1.4 and 1.6, where the initial amplitudes are chosen 
to be Q!i = 6, a2 = 3. The patterns given by fifty random trajectories out of the total 
number of 10,000 analyzed trajectories are shown in Fig. 1.5(b) in the time interval 
gt e (0, 0.5) and Fig. 1.5(c) in gt G (0, 5). The final snapshot of the "cloud" ring at 
£ft — 5 is given in Fig. 1.5(d). 

The method of classical trajectories can be used not only numerically (Figs. 1.5 
and 1.6) but also analytically in special cases. For example, the evolution of low- 
noise fields in the no-energy-transfer regime can be found analytically in the first 
approximation with the solution given by elementary trigonometric functions. To 
show this, let us analyze integrals of motions. On assuming the initial amplitudes of 
the form ai ~ 2r + xi + iyi and 0:2 = r + X2 + iy2, the integrals of motion can be 
expressed in the form of successive corrections 

E = + AEi + AEo^ 

r = 4r3 + AFa + AFi + AFo, (1.26) 

where we denote AEi ~ 4r {xi +2:2), AEq = a;f + + 2 + y|) , and AT2 = 
4r2 (xi + X2) , AFi =r(xl-yl+ 4x1X2) and AFq = X2 {xj - yf) + 2xiyiy2. 
By substituting n2 = E/6 + e, where e is a small correction, and after omitting the 
cubic term 2e^, we find that the denominator in Eq. (1.16) can be approximated by 
the quadratic function 

n2{E -2n2f -V^ !v2E{A'^ -e^) . (1.27) 

Now, we can perform integration of these elementary functions leading to the simple 
result 

n2it) = r'^ + B + Acosngt, (1.28) 

where = V8E,A= {xi - 2x2^ + 3{yl+ yl) and B = |r -h a;2). We 
get a similar result 

ni {t) = E- 2n2 = 4r^ + 4S - 2Acosngt (1.29) 

for the fundamental (or subharmonic) mode. Both solutions are constant functions 
weakly perturbed by harmonic function. The evolution in phase space can be un- 
derstood clearly by analyzing Figs. 1.5 and 1.6. Due to the frequency dispersion 
{{xk,yk}) [see Fig. 1.5(b),(c)], different trajectories are drifting variously and 
create a crescent-shape cloud in phase space, which develops later into a full ring as 
seen in Figs. 1.5 and 1.6. One has to perform the averaging of solutions to calculate 
the required statistical moments. We find ni = 4r^ ,n2 = , and 

^ = 16r'^ + 16W+2'^, 

^ = r^+W+^A^, (1.30) 



12 



gt = gt = 0.5 gt=l 



a 




-5 5 -5 5 -5 5 



gt=1.5 gt = 2 gt = 2.5 




-5 5 -5 5 -5 5 

Re a Re a Re a 



Fig. 1.6 Classical trajectory simulation of quantum evolution of the Q-function for the same 
initial conditions and interaction times as in Fig. 1.4. In our simulation 10,000 trajectories 
were calculated. 

where B = 0,A^ = fr^a^ = ^r^, W = ^r^a^ = fr^ and cos2 Qgt = i. 
Finally, we arrive at the semiclassical Fano factors given by simple rational numbers: 

Ff = 1(41^+11^) =2, 

F| = -L(s^+iA^).^. (131, 

By analyzing Figs. 1.1(b) and 1.3 as well as tables 1.1 and 1.2, we conclude that 
our estimations (1 .3 1) are in very good agreement with those values of Fano factors 
obtained by the quantum numerical analysis of Sect. 1.2.2. 
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1.3 fflGHER-HARMONIC GENERATION 
1.3.1 Quantum analysis 

In this section, we will generalize our results of Sect. 1.2 to describe the processes 
of the A'^-th harmonic generation. Again, we will focus on predictions of the sub- 
Poissonian photon-number statistics. 

Processes of the 7V-th harmonic or subharmonic generation can be described by 
the conventional interaction Hamiltonian (e.g., [5]) 

H = hg (a^ a^'j^ + aY^ ciN^ (1.32) 

for A'^ = 2, 3, • • •. In (1.32), oi and a^v denote annihilation operators of the fun- 
damental and iV-th harmonic modes, respectively, and 5 is a nonlinear coupling 
parameter. For short evolution times, the following approximation of the quantum 
Fano factors can be obtained for the fundamental mode [13] 

= 1-2N{N -1) r^-'^r2 sine gt + 0{{gtf} (1.33) 

with N = 2,3, and for higher harmonics: 

F^ = 1 - 36rfr2 {rj + 2) sm 9 {gtf + 0{{gtf}, 

F^ = 1 -64rir2 {17 + Urf + 2^) sine {gtf + 0{{gtf}, (1.34) 

where rk are input amplitudes, and e = A'^^i — <j)N is the input phase mismatch. 
For spontaneous harmonic generation (i.e., for r^v = 0), Eqs. (1.33)-(1.34) simplify 
to the formulas derived by Kozierowski and Kielich [14]. This analysis shows the 
possibility of sub-Poissonian light generation in short-time regime under the proper 
phase condition. 

On testing different coherent input amplitudes and phases in order to minimize the 
Fano factor for long-interaction times, we have discovered a regime, for which the 
harmonic field exhibits the quasi- stationary sub-Poissonian photocount noise. The 
regime occurs if the ratio of amplitudes \ai\ and |ajv| is equal to N, and phases 
are related by iV^i = 4>n- As described in Sect. 1.2 for SHG, this is a quantum 
analog of the no-energy-transfer regime [25] known from classical nonlinear optics 
as an evolution exhibiting the no-energy transfer between the interacting modes. The 
intensities of both modes remain quasi- stationary during the interaction. Obviously, 
in quantum analysis some small energy fluctuations between modes are observed as 
a consequence of vacuum fluctuations. However, the influence of energy fluctuations 
can be neglected for strong fields. 

For better comparison of theoretical predictions for different order processes, we 
have plotted the quantum Fano factors for both interacting modes in the no-energy- 
transfer regime with iV = 2 — 5 and r = 5 in Fig. 1.7. One can see that all curves 
start from F^j^ (0) = 1 for the input coherent fields and become quasi-stationary 



14 



3.0 



1- 25 

o '^■^ 

-I— • 

o 
ca 

2.0 

O 

c 
ca 



1 .5 



1.0 





(a) 


1.0 
0.9 






0.8 












^ ^ 


0.7 




V 










0.6 








0.5 



10 20 _ 30 



40 



50 




Fig. 1.7 Time evolution of tlie exact quantum Fano factors: (a) = (iV) for tlie 
fundamental mode and (b) F^ for the harmonic mode in ATth-harmonic generation for N=2 
(thickest curve), 3, 4, and 5 (thinnest curve). Time t is rescaled with frequency f2, given by 
(1.52), and coupling constant g. The harmonic-mode amplitude is r = rjv = 5. The dotted 
lines correspond to the semiclassical Fano factors, given by (1.54) and (1.55). It is seen that 
the fundamental mode is super-Poissonian, whereas the harmonic mode is sub-Poissonian for 
all non-zero evolution times. 



after some relaxations. The quantum and semiclassical Fano factors coincide for 
high-intensity fields and longer times, specifically for t > 50/{Clg), where O will be 
defined later by Eq. (1.54). InFig. 1 .7, we observe that all fundamental modes remain 
super-Poissonian (F^ {t) > 1), whereas the A^th harmonics become sub-Poissonian 
(F^ (t) < 1). The most suppressed noise is observed for the third harmonic with 
the Fano factor F^ « 0.81. In Fig. 1.7, we have included the predictions of the 
classical trajectory method (plotted by dotted fines) to show that they properly fit 
the exact quantum results (full curves) for the evolution times t > 50 /{Qg). The 
small residual differences result from the fact that the amplitude r was chosen to be 
relatively small (r = 5). This value does not precisely fulfill the condition r :» 1. 
We have taken r = 5 as a compromise between the asymptotic value r ^ oo 
and computational complexity to manipulate the matrices of dimensions 1000 x 
1000. Unfortunately, we cannot increase amplitude r arbitrary due to computational 
limitations. Numerical values of the quantum Fano factors in comparison with their 
semiclassical approximations for the fundamental mode, given by Eq. (1.56), are 
presented in their dependence on in table 1.1 and Fig. 1.8(a). Analogously, those 
values for harmonics are presented in Fig. 1.8(b) and table 1.2 as calculated by the 
numerical quantum method and from analytical semiclassical formula (1.57). It is 
seen that the approximate predictions of the Fano factors, according to (1.56) and 
(1.57), fit very well the values obtained by applying the numerical quantum method. 
In fact, the differences between the approximate and exact values are hardly visible 
on the scale of Fig. 1.8. Nevertheless, some small (< 1%) differences in F^j^f (see 
tables 1.1 and 1.2) can be explained by the fact that the value of r for numerical 
analysis was chosen too small. 
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1.3.2 Classical analysis 

Our classical analysis of higher-harmonic generation follows the same method as 
described in Sect. 1.2.3. The classical model of the TVth-harmonic generation can be 
described by 

n = g [of + {aD'^UN] , (135) 

which in a special case of = 2 goes over into Eq. (1.11). In Eq. (1.35), ai (a^v) 
is the complex amplitude of the fundamental (A^th-harmonic) mode. Hamiltonian 
(1.35) leads to the pair of complex differential equations [26] 

oiN = -igcti, (1.36) 

On introducing real amplitudes and phases, = rfee"^'', (1.36) can be transformed 
into the system of three real equations: 

ri = —gNr'^^^rN sm9, 
Tn = gvi sin e, 



Table 1.1 Quasi-stationary values of the quantum Fano factors and their semiclas- 
sical approximations , given by (1.54), for the fundamental mode in A^th-harmonic 

generation with N = 1 — 5 in the no-energy-transfer regime. The values of are 
calculated for r = rjv = 5. 

N Ff {F^-Fl)/F^ 

111 

2 1.5029291 3/2 0.0020 

3 1.8202032 29/16 0.0042 

4 2.0323293 101/50 0.0061 

5 2.1830414 13/6 0.0075 



Table 1.2 Same as in table 1, but for the A'^th harmonic mode; F^ are calculated from 
(1.55). 



N 


-''at 


F'^ 


\fS-f^\/fS 


1 


1 


1 





2 


0.83228800 


5/6 


0.0013 


3 


0.81125970 


13/16 


0.0015 


4 


0.81924902 


41/50 


0.00092 


5 


0.83331127 


5/6 


0.000026 
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N N 

Fig. 1.8 Semiclassical (solid bars) and quantum (dithered bars) Fano factors versus order 
N of harmonic generation for (a) fundamental and (b) A'^th-harmonic modes in the quasi- 
stationary no-energy-transfcr regime. Figures (a) and (b) for = 1 — 5 correspond to tables 
1 and 2, respectively. It is seen that the quantum results are well fitted by the semiclassical 
Fano factors. According to both analyses, the third-harmonic mode has the most suppressed 
photocount noise. 

= 5 (rf/rjv -A^Vf-V) cos ^, (1.37) 

where 6 = N(j}\ — (p^ is the phase mismatch. Equations (1.37) have two integrals 
of motion: 

E = r1 + Nr% = m + Nhn, 

r = rfrjvcos6». (1.38) 

On extraction of ri and 6 from Eq. (1.37), we find equation for the amplitude r^: 

{rNTN/gf + =rl{E-Nrl)'' (1.39) 

or its simpler form for the intensity Un = rlf-. 

{hN/2gf = nN{E-NnNf -V". (1.40) 

The general solution for tin [t] is a periodic function oscillating between the values 
Umin and nmax- The solution can be given in terms of the Jacobi elliptic functions 
for N = 2 and N = 3, and in terms of hyperelUptic functions for N > 3. 

One elementary solution of set of Eqs. (1.37) is obtained for the zero initial phase 
mismatch 9 = and the initial ampUtudes satisfying the condition n = Ntn- The 
solution reads as 

ai{t) = ricxp(-i.gtrf"^), 

a^it) = rNCxp{-iNgtr^-^) , (1.41) 

which corresponds to the no-energy-transfer regime, since the amplitude and energy 
in both the interacting modes remain constant m (t) = |q!i (t)| = and riAr (t) = 
|ajv (i)l' = r% [25]. 
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1.3.3 Classical trajectory analysis 



The results of Sect. 1 .3.2 can be used in the method of classical trajectories in analogy 
with the technique described in Sect. 1.2.3. We need to express the trajectories in 
their dependence on small noise parameters Xk and yk- The integrals of motion, given 
by (1.38), can be expressed in a form of corrections in successive powers of large r: 



E = N{N + l)r^ + AEi + AEq, 



where 



AEi 
AEo 



2N (xi +XN)r, 

xl + yl + N {x% + y%) 



and 
where 



r = 7V^r^+i + ATn + ATn-i + ArN-2 + • • - 



(1.42) 

(1.43) 
(1.44) 



ATn 

Arjv-i 



{xi + xn) (Nr) 



N 



N-1 



{^i -yl) + N {xiXN + yiyn) 



{Nrf-\ (1.45) 



The lower-order terms ATis[_2, AFjv-s, ■•• can be neglected in further considera- 
tions. On assumption of high-intensity fields (r ^ 1), we can substitute 



riN 



E 



N{N+1) 



(1.46) 



where e is a small correction of stationary value. Then, r.h.s. of (1.40) can be 
rewritten as 



njv {E - Nhn)^ - r 



N 



N 



\JV-2 



E 



N-1 



(1.47) 



2(iV-|- 1)' 

on omission of higher-order terms involving e^, e^, • • •. One arrives at simple equation 



■,N-2 



E 



JV-l 



2{N+1) 
Thus, the solution of (1.40) reads as 

n;v (t) =r'^ + B + AsinClgt, 



(1.48) 



(1.49) 
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where the frequency is given by 



and 



A = 



n = J ^ (1.50) 



-\/4{xi - NxNf + 2NiN + 1) (yi - yjv)' 



TV- 

From (1.38), a result similar to (1.49) is obtained for the fundamental mode: 

ni {t)=E- Nun (t) = N^r^ + N^B - NAsinngt. (1.52) 

It is seen that both solutions (1.49) and (1.52) are given in a form of large constants 
weakly perturbed by harmonic function. 

Now, on applying the classical trajectory method, one should perform averaging 
over all solutions (1.49) and (1.52) to calculate the required statistical moments. 
Here, we calculate the first and second-order field-intensity moments necessary for 
determination of the Fano factors. The mean intensities of the fundamental and 
harmonic modes are simply given by ni = N'^r'^ and un = r^, respectively. The 
second-order moments of field intensity are found to be 

^ = N'^r^ + N'^W+^N'^A^, 

^ = r^ + s2+l]42. (153) 

in terms ofA^ = r^ {2 N'^ + N+ 1)/{N +lf andW = 2r^ /{N +lf . We note that 
B vanishes. The term sin^ V.gt can simply be estimated with 1 /2 for sufficiently long 
time t, when Uk (t) and (t) become quasi- stationary. The relaxation in rife (t) and 
(t) is observed for short times t due to the presence of harmonic sine function 
and residual phase synchronization. The mean value of the frequency (1.50), given 

by _ 

O « ^/2N {N + 1) (ATr)^-^ , (1.54) 

enables estimation of the oscillation period Tqsc = 2-k/Q., whereas the standard 
deviation 

Al^ « v'2iV(iV + l)iV^-V^-2^^ (1.55) 

determines the duration Troi = 27r/Af2 of relaxation. By comparing the characteristic 
times Tqsc and T^ei, one finds that the evolution time can be scaled by r = ilgt to 
synchronize optimally the oscillations of the exact quantum solutions for different 
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N. These synchronized oscillations of the Fano factors are clearly presented in Fig. 
1.7. 

Finally, we arrive at the semiclassical Fano factors 



which are the compact-form analogs of the quantum Fano factors. The semiclassical 
Fano factors for the fundamental and higher harmonics for various values of N are 
Usted in tables 1.1 and 1.2, and plotted in Figs. 1.7(a) and 1.7(b), respectively. 

Oursolutions(1.56)and(1.57)reducetotheresultsderivedinRef. [19]forA^ = 2, 
and those of Ref. [20] for = 3. By analyzing (1.57), we find that higher harmonics 
evolve into quasi- stationary sub-Poissonian states (F^ < 1) for any N > 1. Except 
for second harmonic, the photocount noise reduction in higher harmonics becomes 
less effective with increasing N. Thus, the deepest noise reduction occurs for the third 
harmonic as described by the Fano factor Ff = j| =0.8125. The photocount noise 
reductions for the second and fifth harmonics are predicted to be the same, although 
the quantum analysis (see table 1.2) reveals that they differ slightly (< 1%). As 
comes from (1 .56), the fundamental mode has solely the super-Poissonian photocount 
statistics (Ff > 1) with noise monotonically growing in N for the no-energy-transfer 
regime. For N = 1, the process is linear and no change in the photon statistics occurs. 
The interacting modes remain coherent with the unit Fano factors for both modes. It 
is worth noting that qualitatively different photocount statistics of the fundamental 
mode is observed in the short-interaction regime as given by Eqs. (1.6) and (1.33)- 
(1.34). 

We have shown, in agreement with the results presented in [ 1 9] , that the method of 
classical trajectories gives very good predictions in the case of strong-field interactions 
(i.e., for the photon numbers larger than 10). The calculation speed of the method 
does not depend on numbers of interacting photons. But better approximation is 
achieved with the increasing number of photons. Thus, the method is very fast and 
significantly simplifies the tedious exact quantum calculations. 

1.4 CONCLUSION 

We have presented quantum, classical and semiclassical descriptions of second and 
higher harmonic generations. We have demonstrated that these processes can be a 
source of sub-Poissonian light. On testing different coherent input amplitudes and 
phases in order to minimize the Fano factor, we have discovered a quantum regime 
for which the long-interaction output is generated with the quasi- stationary sub- 
Poissonian photocount noise [19, 20, 21]. The regime occurs if the initial coherent 
state amplitudes are related by |ai| = A''|Q!jv| andArg(Q!jv) — A''Arg(ai). Thisisa 



1 6N^ + N + 1 

2 {N + lf ' 

1 27V2 + N + 5 

2 (N + lf ' 



(1.56) 
(1.57) 
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quantum analog of the no-energy-transfer regime [25] known in classical nonlinear 
optics as an evolution exhibiting no-energy transfer between the interacting modes. 
The intensities of both modes remain quasi-constant in time during the interaction. 
Obviously, in a quantum analysis some small energy fluctuations between modes are 
observed as a consequence of vacuum fluctuations. However, the influence of energy 
fluctuations can be neglected for strong fields. 

We have proved that in the no-energy-transfer regime, the fundamental mode 
evolves into a quasi- stationary state with the super-Poissonian {F^ > 1) photocount 
statistics, whereas the A'"-th harmonic goes over into a sub-Poissonian {F^ < 1) 
quasi-stationary state. We have found that the most suppressed photocount noise 
is obtained for the third harmonic as described by the quantum Fano factor = 
0.811 • • •. Good analytical predictions of the quantum Fano factors for both the 
fundamental and harmonic modes (Fg = 13/16 = 0.8125) were obtained under the 
semiclassical approximation in the strong-field limit. 
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